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The combination of spin transitions and magnetic ordering provides an interesting struc- 
ture of phase transitions in Prussian blue analogues (PBAs). To understand the structure of 
stable and metastable states of Co-Fe PBA, it is necessary to clarify free energy as a function 
of magnetization and the fraction of the high-temperature component. Including the mag- 
netic interaction between high-temperature states, we study the magnetic phase transition 
of Co-Fe PBA in addition to spin transitions. Here, we take into account the degeneracy 
changes due to charge transfer between Co and Fe atoms accompanying the spin transition. 
In this study, the charge transfer between Co and Fe atoms is explicitly taken into account 
and also the ferrimagnetic structure of Co-Fe PBAs is expressed in the proper way. First, 
we found systematic changes in the structures of stable and metastable states as functions 
of system parameters using mean field theory. In particular, the existence of a metastable 
magnetic-ordered high-temperature state is confirmed at temperatures lower than that of 
the hysteresis region of spin transitions. Second, we found that the magnetic interaction 
causes complex ordering processes of a spin transition and a magnetic phase transition. The 
effect of a magnetic field on the phase structure is also investigated and we found metamag- 
netic magnetization processes. Finally, the dynamical properties of this metastable state are 
studied by Monte Carlo method. 

KEYWORDS: spin transition, charge transfer system, metastable structure, ferrimagnetism, 
metamagnetism 

1. Introduction 

In some octahedral coordinate iron complexes, a spin crossover (SC) may occur between 
a high-spin (HS) state and a low-spin (LS) state. An SC transition is induced by external 
stimulations such as temperature, photoirradiation, and mag netic field. 1 - 19 ' 21 - 31 ) The HS state 
is favorable at high temperature owing to its high degeneracy. On the other hand, the LS state 
is favorable at low temperature because it has a low energy with a low degeneracy. Interesting 



*E-mail address: konismQspin.phys.s.u-tokyo. ac.jp 



1/17 



J. Phys. Soc. Jpn. Full Paper 

transitions between HS and LS states are provided by the competition between entropy gain 
and energy gain. It has been pointed out that cooperative interactions are important for SC 
transitions. According to the system parameters, smooth transition or discontinuous first- 
order phase transition occurs in the SC transition. 4 ~ 4 ' 15 ~ 23 ) This situation is well described 
by the Wajnflasz-Pick (WP) model, 15 ) which gave the theoretical basis of the mechanism 
of the cooperative transitions using an Ising model with degenerate states. The WP model 
and its extended models have explained successfully the static and dynamical properties of 
various cases of SC transitions. 19 ~ 23 ) Control between the HS and LS states has been realized 
by photoirradiation as light-induced excited spin state trapping (TIESST), 7 " 12 ' 19 - 24 ) and the 
structure of the photoinduced metastable state of the systems has become an important 
topic. 11 - 14 ) 

Recently, we have reported that the metastable state can exist intrinsically at low tem- 
perature, which has been supported experimentally. 13 ' 14 ) The effects of applying a magnetic 
field have been studied in some SC complexes. 25 ~ 31 ) Since the magnetic moment in the HS 
state is larger than that in the LS state, the HS state is stabilized by applying a magnetic 
field. A shift in transition temperature under an applied magnetic field was confirmed exper- 
imentally. 25 ~ 27 ) Using a pulsed high magnetic field, the creation of the HS state from the LS 
branch in the thermal hysteresis loop was observed and these results were interpreted using 
an Ising-like model. 28 ~ 31 ) However, since both the HS and LS states in the SC complexes are 
paramagnetic, the study of the effects of applying a magnetic field has been performed only 
for the paramagnetic region. 

Some of the Prussian blue analogues (PBAs) 32 ~ 48 ) show thermal-induced phase transi- 
tions and photomagnetic effects. 39-45 ) In such materials, a strong cooperativity operates since 
transition-metal ions are bridged by CN ligands with a three-dimensional network structure. 
Co-Fe PBA is an attractive material owing to its two-way photoswitching, i.e., magnetic <-> 
nonmagnetic. 39 ~ 42 ) The thermal induced phase transition between two states, i.e., Co(II)(HS, 
S = 3/2)-Fe(III)(LS, S = 1/2) and Co(III)(LS, S = 0)-Fe(II)(LS, S = 0) is called a charge- 
transfer-induced spin transition (CTIST) 42 ) in Co-Fe PBA. Since the state of Co(II)-Fe(III) 
is favorable at high temperature owing to its high degeneracy, we call it a high-temperature 
(HT) state. Since the Co(III)-Fe(II) state is favorable at low temperature owing to its low en- 
ergy, we call it a low-temperature (LT) state. In this material, the magnetic coupling between 
Co and Fe ions is antiferromagnetic. Therefore, the HT state is ferrimagnetic. 

In this work, we study the effects of magnetic interaction and magnetic field on the co- 
operative phenomena of spin-crossover transitions and magnetic ordering for a Co-Fe PBA 
system using an extended Wajnflasz model. The calculations are performed using the Mean 
Field (MF) approximation and Monte Carlo (MC) methods. In particular, we study the effects 
for the metastable HT state ,which exists in low temperature region. 
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The structure of this paper is as follows. In § 2, the model used in this study is explained. 
In § 3, the MF approximation used in this study is explained. In § 4, the effects of magnetic 
interaction and applying a magnetic field are presented using the MF approximation. The 
calculation results using the MC method are presented in § 5. Section 6 is devoted to Summary 
and Discussion. 




Fig. 1. Schematic illustration of present model of Co-Fe PBA. Large blue circles are Co ions and red 
circles are Fe ions. Small white circles are N atoms and gray circles are C atoms. 



2. Model 

Let us consider the CTIST phenomenon of Co-Fe PBA. This material consists of a bipartite 
lattice. One of the sublattices is occupied by Fe ions and the other is occupied by Co ions 
(Fig. 1). We call a site of the Fe sublattice an A-site, and a site of the Co sublattice a B-site. 
The LT state consists mainly of Fe(II) and Co(III) ions. The HT state consists mainly of 
Fe(III) and Co(II) ions. Here, we consider the degeneracy of the spin degree of freedom. At an 
A-site, Fe(II) in the LT state is S = 0, whose degeneracy is 1, and Fe(III) in the HT state is 
5=1/2 with a degeneracy of 2. At a B-site, the degeneracy of n-Q = is 1 and that of n-Q = 1 
is 4 because Co(III) in the LT state is S = 0, and Co(II) in the HT state is S = 3/2 with a 
degeneracy of 4. These degeneracies due to the spin states are listed in Table I. Besides the 
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g p honon(LS). The degeneracy due to phonons is much larger than that of spin. However, in 
this study, we mainly use spin degeneracy, which is enough to provide the general aspects of 
the phase structure. 

This charge transfer phenomenon is expressed by an electron transfer between Fe and Co 
atoms. To express this transfer, we introduce the quantity n, which is 1 for Fe in the LT 
state (i.e., riA = 1) and in the HT state (tia = 0). Correspondingly, it is for Co in the 
LT state (ne = 0) and 1 in the HT state (ne = 1)- Here, denotes A-(B-)sites. In this 
system, electrons are transferred between A and B sites. We have introduced the following 
Hamiltonian for this system: 

H = Dj2ni + eJ2nini + 4jJ2s\si-2h I J>A + E 4 I ■ « 

3 \ i 3 J 

Here, D (> 0) is the difference in on-site energy between A- and B-sites, e represents the 
repulsive interaction between the electrons at nearest-neighbor sites, and J represents the 
magnetic interaction between A- and B-sites, and s denotes the spin. The external magnetic 
field is given by h. Since an Fe and Co pair has an electron, the following relation holds: 

E< + E n B = ^ (2) 

» 3 

where N is half of the total number of Co and Fe ions of the lattice. Here, we also introduce 
the magnetic states of the spins by sa and sb- When ha = 1, S = and thus sa takes only 0, 
and when ha = 0, S = 1/2 and thus sa takes —1/2 or 1/2. Similarly, when hb = 0, S = and 
sb takes 0, and when hb = 1, S = 3/2 and sb takes —3/2, —1/2, 1/2 or 1/2. Specifying sa and 
sb, the degeneracy is taken into account naturally. In this model, electrons tend to stay at 
A-sites energetically, which denotes the LT state. When all the electrons are at A-sites, there 
is no energy cost due to D. We define this state as the perfect LT state. At high temperature, 
electrons tend to stay at B-sites owing to the difference in degeneracy mentioned above, which 
corresponds to the HT state. 

Let us briefly review the properties of this type of Hamiltonian in the case with neither 
magnetic interactions nor a magnetic field (J = h = 0). 13 ) In this case, it is known that the 
critical temperature (7i/ 2 ) of charge order, at which (hb) = 1/2, is exactly given by 

kBT 1/2 = ^15^ (3) 

where /cb is the Boltzmann constant. <7ht and <7lt are the degeneracies of the HT (= 8) and 
LT (= 1) states, respectively. With these <7ht and ^lt, we have 

k B T 1/2 = D/ln8 ~ 0.4809D. (4) 

At this temperature, (wb) changes discontinuously in the first-order phase transition, when 
T1/2 is larger than a critical value Ti s i ng /4, where 7i s i ng is the critical temperature of the Ising 
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model on the same lattice. On the other hand, if Ti s i ng /4 < Ti/ 2 , (^b) changes smoothly. On 
the simple cubic lattice, /cB^ising = 4.51e. Using this value, whether the transition is smooth 
or discontinuous depends on the ratio e/D, and the critical ratio is 

— = ~ 0.4266. (5) 

D 4.51 In 8 V ' 

In the mean field theory, Ti s i ng MF = 6e, and we have 

€r 4 



D 61n8^°- 3206 ' (6) 
which is different from eq. (5). Furthermore it has been pointed out that there are other 
critical values of e. When e exceeds a critical value 

e CG /D = 1/6, (7) 

a low-temperature metastable branch appears. Moreover, when e is larger than a critical value 



ecx _ 1 1 + \J 9ht/ 9lt ^ 
D 6 1 - v / 3ht/ / S'lt' 
a local stability of the HT state remains down to T = 0, i.e., the low-temperature limit of the 
hysteresis loop reaches T = 0. 

3. Mean Field Approximation 

In this section, we present a mean field formulation of the system with external magnetic 
interaction (J) and the applied magnetic field (h). To study the model in the grand canonical 
ensemble, we extend the Hamiltonian including the chemical potential fx. 

H = D^n B +e^ n\n 3 B +4J £ s\s^-2h I £ s\ + £ 4 j ( £ n\ + £ n{ j . (9) 

In the mean field theory, the Hamiltonian is approximated as 

Hmft = (6n B e - fj) n A + [6«A£ + £> - ^] ^ n B + (24Js A - 2/i) ^ s B 

i i i 

+(24Js B - 2h) ^ A - 6Nen B n A - 24NJs A s B , (10) 

i 

where n A , n B , s A and s B are the averages of quantities described as 

The partition function is given as 

Z = Trexp(-/3W M FT) 

= ( e -/*(6™-f) + 2 cosh (/3(12Js B - h))) N 
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x (1 + 2e- (3l6nAt+D -^ [cosh (/3(12 Js A - h)) + cosh (3/3(12 Js A - h))]) N . (12) 

Here, it should be noted that Tr denotes the summation over all states, (n A ,s A ) = 
(1, 0), (0, -1/2), (0, 1/2) and (n B , 4) = (0, 0), (1, -3/2), (1, -1/2), (1, 1/2), (1, 3/2). 
Using the Hamiltonian (10), we have 



2e- /3 ( 6nA£+D -^[cosh(/3(12Js A - h)) + cosh(3/3(12Js A - h))] 
(n B ) - n B - — 2e -/3(6n Ae +D- A1 )[ cosh ( j g( 12 j SA _ ft)) + cos h(3/?(12JsA - fc))] ' 

and 



e -/3(6n B e-^) 

<n A ) = n A = ^7^Z7^ - 2 cosh(/3(12 j SB _ h)) ■ ( 14 ) 

Here, because of the condition for the electron number in eq. (2), n A + n-g = 1. Hereafter, we 
express the HT fraction 

n = ne, (15) 

and thus 

n A = 1 — n. (16) 

Then, we have 

2e -/3[6(l-n)e+D- l x] [ cos h(/3(12 Js A - /l)) + COsll(3/?(12 Js A - /l))] 

1 + 2e-^( 1 - n ')"+ D -^ [cosh(/3(12 Js A - /i)) + cosh(3/3(12Js A - h))} 

e -/3(6ne-^) 

+ e _^ (6ne _^) - 2cosh(/3(12Js B - /*)) ~ ^ 
from which we determine the chemical potential \i. Setting / = e 13 ^, we obtain 



f = J cosh(/3(12Js B - h)) «d +3e) 

J y cosh(/3(12Js A - ft)) + cosh(3/3(12Js A - h)) ' y 1 

The free energy F is given by 

F = -llaZ 

= -Tln(/e~ 6/3en + 2cosh(/3(12Js B -/»))) 

-Tin [1 + 2/e^ /3 ( 6e ( 1 - n ) +D )(cosh(/3(12Js A - /i)) + cosh(3/3(12Js A - h)))} 

-6en(l - n) - 24Js A s B - (19) 

The self-consistent equations for s\ and S2 are obtained by the conditions of dF/ds A = 
dF/dsB = 0. The condition OF/Osa = leads to 

sinh(/3(12J.s A - h)) + 3sinh(3/?(12Js A - h)) 
SB ~ ~ cosh(/3(12Js A - h)) + cosh(3/3(12Js A - h)) ( ^ 

and from OF/Osa = 0, we have 

sa = - tanh(/?(12 Js B - h))n. (21) 
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The magnetization of the system is given by 

m = SA + SB- (22) 

To obtain the thermal properties of the system, we have to find solutions that satisfy the 
eq. (13), (20) and (21). 

4. Phase Transition in Mean Field Theory 

Hereafter, we take D as the unit of energy. First, we review the properties of the system 
without magnetic interaction (D,J,h) = (1,0,0) as a reference. Let us study the effect of 
repulsive interaction, e. The temperature dependences of the HT fraction (n(T)) obtained by 
solving eq. (13) with eq. (18) are presented in Fig. 2(a) for the values of e = 0,0.1,0.2,0.34, 
and 0.4. When e = 0, the transition between HT and LT states is gradual and continuous. 
Until the critical value of e(= 0.32D), eq. (6), the transition is continuous. When e > e c , the 
transition becomes discontinuous and it shows a first-order phase transition. In this case, the 
transition temperature (T c ) is T x i 2 = 0.48 (eq. (4)), and the system has three solutions at T c . 
The solutions with the largest and smallest HT fractions are stable states, whereas that with 
the intermediate HT fraction represents an unstable state. For ecc = 1/6 < e < ecx — 0.349, 
the HT state is the metastable state, as shown in Fig. 2(a). For large values of e (e > ecx), 
the local minimum of the HT state exists at all temperatures as shown in Fig. 2(a) (e = 0.4). 
In Fig. 2(b) the corresponding free energy profiles for e = 0.34 are depicted for various 
temperatures, where we find the metastable structure at low temperature. 
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Fig. 2. (a) Temperature dependence of HT fraction with parameters of J = 0, D = 1 and e = 
0, 0.1, 0.2, 0.34 and 0.4. (b) Free energy curves as functions of HT fraction. 
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4-1 Effect of magnetic interaction 

Let us consider the effect of magnetic interaction using J without the external magnetic 
field {h = 0). The temperature dependences of n are depicted in Fig. 3(a) with the parameters 
(D,e,h) = (1, 0.34, 0) with J = 0, 0.02, 0.04, and 0.06, where the system possesses a HT 
metastable branch at low temperature in the case without the magnetic interaction (J = 0). 
When the magnetic interaction is included, we find that the metastable branch is enlarged. 
In the case of J = 0.02, the HT fraction of the metastable HT state is increased. The HT 
fraction of the unstable solution is also affected. Here, it should be noted that the solution 
of J = is not changed unless the magnetization has nonzero solution in the MF treatment. 
Indeed, the thermal hysteresis loop does not change at all in this case. When J is increased 
to 0.04, the metastable branch is further enlarged and terminates at T = 0.41, at which 
n changes discontinuously to the solution of J = (black line), where the magnetization 
vanishes discontinuously. We call this temperature T cmag . The temperature dependences of 
magnetization m(T) are depicted in Fig. 3(b). For J = 0.06, the metastable branch terminates 
at T = 0.57. The magnetic states in the case J = 0.02 and J = 0.04 are metastable, whereas 
it is stable, i.e., the equilibrium state in the case J = 0.06. This means that the model is a 
simple ferrimagnetic without CTIST. 

In Fig. 3(a), we find that the metastable solution (solid green line) and the unstable 
solution (dotted green line) terminate at different temperatures T = 0.41 and T = 0.39, 
respectively, and are not connected to each other. This is due to the fact that the present 
model has three order parameters, i.e., (n, sa, sb)- At T = 0.35, there are three solutions that 
are shown in Fig. 4(a) by closed circles for local minimum points and an open circle for a 
saddle point. When T is increased up to T = 0.4, the saddle point disappears, as shown in 
Fig. 4(b), where two stable solutions correspond to points on the green solid line and black 
solid line denoting the LT state. At T = 0.42, the metastable point disappear and only one 
stable point remains, which gives a point on the black line. For large values of J, HT remains 
locally stable until T = 0, as we see in Fig. 3(a) for J = 0.06. 

Here, we point out a unique dependence of the stability of the ferrimagnetic state on 
J. If J < J c = D/18, the ferrimagnetic state is metastable and not an equilibrium state 
at low temperatures. Here, J c is obtained by the comparison of the energy of the complete 
ferrimagnetic state and the LT state. In Fig. 5, we show an example of J = 0.54 < J c . There, 
the equilibrium solution is drawn by a bold solid curve. 

4-2 Effect of external field 

Let us consider the effect of the applied magnetic field h. Here, we study the case of 
e = 0.34 where both thermal hysteresis and the HT metastable exist at low temperature. The 
magnetic field dependence of magnetization at T = 0.1 is presented in Fig. 6(a) for various 
values of J (J = 0.01,0.04, and 0.06), where the T cmag values of J = 0.01,0.04, and 0.06 
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Fig. 3. (a) Temperature dependence of HT fraction for e = 0.34, D = 1 and J = 0, 0.02, 0.04 and 
0.06. (b) Temperature dependences of magnetization for e = 0.34, D — 1, and J = 0.02,0.04 and 
0.06. The solid curves represent stable and metastable solutions and the dashed curves represent 
unstable solutions. 



are 0.133, 0.39, and 0.57, respectively. In the case of J = 0.01, the LT state is a stable state 
in the absence of a magnetic field. For a weak field, the system is always in the LT phase. 
The magnetization is suddenly induced up to nearly the saturated value at h = 0.338. The 
magnetization value of 2 indicates the ferromagnetic state of Fe (S = 1/2) and Co (S = 3/2) 
ions in the HT state. In the case of J = 0.04, a stable state in the absence of a magnetic 
field is still the LT state. At h = 0.281, the magnetization is induced up to 1 and shows a 
plateau until h ~ 0.45. After that, the magnetization is induced up to 2. The magnetization of 
1 indicates an antiferromagnetic coupling between Fe and Co ions in the HT state, indicating 
that the plateau at the magnetization of 1 is the region of the ferrimagnetic state due to the 
antiferromagnetic interaction J. For the case of J = 0.06, the ground state is ferrimagnetic 
even at h = 0. The magnetization is induced up to 2 in a strong magnetic field. These results 
indicate that the applied magnetic field can induce transitions between the three states: the 
LT state, the antiferromagnetic coupling HT state, and the ferromagnetic coupling HT state. 

Next, we study the temperature dependence of magnetization. Figure 6(b) shows the 
magnetization with J = 0.04 at T = 0.1, 0.4, and 0.47. In the case of T = 0.4, the system is 
in the LT state in a small field and the magnetization gradually increases until h = 0.175. In 
this field, the magnetization jumps to 0.87, and then gradually increases up to 2. In this case, 
the metastable state exists, indicating that the hysteresis loop can be observed. In the case 
of T = 0.47, the system is in the LT state in a small field and the magnetization increases up 
to 2 with a jump at h = 0.090. 

Now, we study the effect of the magnetic field on the structure of n(T). In Fig. 7, n(T) 
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Fig. 4. Contour plots of free energy at e = 0.34, D = 1 and J = 0.04. Closed circles indicate the local 
minimum points and the open circle indicates the saddle point. 



are presented for (D,e,J) = (1, 0.34, 0.04) and h = 0,0.01,0.05,0.1, and 0.2. As h increases, 
the HT fraction at high temperature increases. In the metastable HT state, the HT fraction is 
increased. This tendency is enhanced as h increases. In contrast to that in the case of Fig. 3(a) 
where J changes, the HT fraction changes even above T cmag , in the present figure. 

5. Monte Carlo Simulations 

In the previous section, we study the static properties of the model using the mean field 
approximation, where the effect of short-range fluctuation is not taken into account. There- 
fore, in this section, we study the model by Monte Carlo simulation and confirm the results 
obtained in the previous section. Furthermore, we study dynamical properties and discuss the 
metastability of the model. 
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Fig. 5. Temperature dependence of HT fraction for D = 1, e = 0.34 and J = 0.054. The dotted 
curves indicate solutions of the self-consistent equations for J = 0. The solid curve indicates the 
stable solution for J = 0.054. 




Fig. 6. Magnetization curves for (a) T = 0.1, J = 0.01,0.04,0.06, (b) J = 0.04, T = 0.1,0.4,0.47. 
The bold curves indicate the stable solutions and the dashed curves indicate the metastable and 
unstable solutions. 



5.1 Equilibrium states 

Here, we study the changes in n(T) and m(T) starting from the initial states of n = 1 and 
m = 1, respectively, which correspond to the photoinduced HT saturated state at very low 
temperature. At each step, we performed 10000 MCSs for transient steps and 20000 MCSs to 
measure the physical quantities. The system size is 16 x 16 x 16, which is large enough to study 
thermal properties. We heat the temperature up to T = 1 and then cool it down to the original 
temperature. We investigate n(T) in the case of (D, e) = (1, 0.4) and J = 0, 0.025, 0.05. Results 
of these simulations are depicted in Fig. 8. In the case of J = 0, we found a smooth n(T) in 
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Fig. 7. Temperature dependences of HT fraction for D = 1, e = 0.34, J = 0.04 and h = 
0,0.01,0.05,0.1 and 0.2. 



Monte Carlo simulation, whereas, in MF theory, we find that the HT phase exists down to 
T = as a metastable state that causes the first-order phase transition at some temperature 
(Fig. 2(a) e = 0.4). This fact indicates that in MC, different types of n(T) appear from the 
MF result for the same set of parameters, which is naturally expected. Nevertheless, in MC, 
we find a similar trend of n(T) to that of MF (Fig. 3(a)). Here, it should be noted that in MC, 
we find n(T) at high temperature, where the state is paramagnetic, is also caused by J. In 
the MF result, n(T) did not change up to m = 0. Indeed, when J = 0.025, n(T) changes from 
that for J = 0, although the system is paramagnetic in both cases. If we increase J further, 
the system shows a first order phase transition. This trend is the same as that in MF theory. 

Because we sweep the temperature, the result depends on sweeping rate. However, we 
find that the position of the transition from metastable HT to stable LT changes little with 
sweeping rate (Fig. 9), which suggests that the state is stable or metastable and is rather well 
defined by the system parameters. 

5.2 Relaxation and evaluation of unstable states 

Finally, we come to the structure of the metastable state, 13 ) e.g., the local minimum of 
this metastable HT state exists at low temperature. The existence or nonexistence of the 
metastable phase is a significant issue in the field of photoinduced phase transition because 
whether the photoinduced HT phase is a metastable phase or a photogenerated thermodynam- 
ically unstable phase is of current interest. From this point of view, we study the relaxation 
processes from the metastable HT to LT states by MC simulation. Here, we adopt the initial 
state of n — 1 where the positions of HT sites are chosen randomly. Such a state can be 
produced by rapid cooling or light irradiation where the HT sites locate uncorrelatedly. Fig- 
ure 10(a) shows the relaxation processes. In this section, we study the system using (D, e, J, h) 
= (1, 0.4, 0.025, 0) at T = 0.1 for different HT fractions, e.g., 0.87, 0.88, 1.00. The system 
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Fig. 8. MC simulations of temperature dependences of HT fraction and magnetization for D = 1, 
e = 0.4, and J = 0,0.025,0.05. 
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Fig. 9. MC simulations of n(T) for D = 1, e = 0.4 and J = 0.025. We change sweeping rate from 10 4 
MCSs to 10 5 MCSs. The transition temperature is around T = 0.18 for all sweeping rates. 



size is 128 x 128 x 128. When the HT fraction is smaller than 0.90, the relaxation curves show 
monotonic decreases. In contrast, when n is larger than 0.90, the relaxation curves increase 
initially, as shown in Fig. 10(b). This is a particular phenomenon due to the existence of a local 
minimum of the HT state. After that, the states relax to the LT state with sigmoidal shapes. 
Here, we find a separatrix at n sp = 0.89 that denotes the position of an unstable point in MC. 
We study the position of n sp at other temperatures and estimate the temperature depen- 
dences of unstable points, which are plotted in Fig. 11 as dots with a stable n(T) dependence. 
The stable and unstable points show a qualitatively similar temperature dependence to that 
obtained by MF. 
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Fig. 10. Relaxation of HT fraction for D = 1, e = 0.4, J = 0.025 and T — 0.1. (a) These simulations 
are started from n = 0.87,0.88, 1.00. (b) Magnified figure for short time. 
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Fig. 11. Stable and unstable points for D = 1, e = 0.4 and J = 0.025. The squares indicate the 
unstable points determined by MC simulation. The circles indicate equilibrium points. 



6. Summary and Discussion 

The spin states and magnetic properties of Co-Fe PBA were studied. We introduced a 
model in which we can take into account the fact that charge transfer occurs between Co and 
Fe ions that are ordered antiferromagntically at low temperatures. The effects of magnetic 
interaction and an external magnetic field on the structure of the metastable states of n(T) 
and m(T) were investigated by MF and MC methods. We found a systematic change in the 
temperature dependence of n(T) on the parameters of the system, which was found in our 
previous study. Namely, in that study, by changing the parameter e, we reproduced the same 
trend of n(T). In this work, in particular, the combination of the ordering processes of n(T) 
and m(T) is studied. The magnetic interaction induces a magnetic order and enhances the 
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ordering of n(T). The low-temperature metastable HT branch is enlarged and even connected 
to the HT branch. Furthermore, we found a unique temperature dependence of the set of 
(n(T) , m(T)) . That is, the solution of (n(T),m(T)) disappears discontinuously, which causes 
a discontinuous transition of n(T) and m(T). The external magnetic field also causes a change 
in n(T) in a similar sequence to that found for magnetic interaction. In our formulation, the 
ferrimagnetic structure is taken into account, and thus we found that the magnetization shows 
two steps, i.e., among the LT state, the ferrimagnetic state and the field-induced ferromagnetic 
state. The qualitative features of the effect of magnetic interaction obtained using the MF 
approximation were confirmed by MC simulation, although in the MC method we found 
that the magnetic interaction has an effect even in the paramagnetic states. In order to 
characterize the low-temperature metastable branch, we studied the relaxation from the HT 
state. We found that the direction of initial relaxation can be upwards or downwards when the 
metastability exists depending on the initial n. That is, n relaxes monotonically to LT when 
n is smaller than the threshold value n sp , whereas it first increases to its metastable value 
when n > n sp . Studying the temperature dependence of n sp , we reproduced a temperature 
dependence similar to that of the unstable solution in the MF study. 

The magnetic properties of PBA obtained in this study, particularly the trend of the 
changes in the temperature dependences of n and m on system parameters will be useful for 
classifying various materials, not only those belonging to Fe-Co PBA but also more general 
materials in which charge transfer and spin-crossover induce phase transitions. In particular, 
in the present formalism, we can take into account details of microscopic structure changes, 
such as the spin values in the HT structure, and thus the present formalism will be useful for 
studying the combined phenomena of charge transfer, spin-crossover, and magnetic ordering 
in a wide range of materials. 

In the present study, we use a short-range interaction e as an elastic interaction approxi- 
mately. However, to be realistic, this elastic interaction is long-range. In this type of materials, 
this long-range interaction is important to study the effect of volume change or that of pres- 
sure. The effect of long-range elastic interaction will be reported elsewhere in the near future. 
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